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Abstract 

We describe a general methods to localize any sort of k-separability and therefore also the corresponding 
partial entanglement in genuinely multipartite mixed quantum states. Our methods are based exclusively 
on the known twopartite methods and some systematic procedures of combinatorial nature. Our methods 
are formalized in an algorithmic-like fashion and therefore they are easily implementable in a computer 
^ • environments and might be effectively used for studying numerically these questions for sufBciently small 

systems 

^ ■ 1 Introduction 

Quantum entanglement, since its discovery by Einstein, Podolsky and Rosen [T], and Schrodinger [2] at- 
tracted very large attention, especially in the past two decades due to its intriguing properties and extremely 
important significance in the quantum information processing tasks [3], [4], [11) . The main efforts have been 
focused on bipartite entanglement and at present many interesting and deep results have been obtained in 
this case, see [5], [6] for some recent reviews. 

The present note is addressed to a very preliminary discussion of certain aspects of multipartite entan- 
glement. The importance of multipartite entanglement can be illustrated by several illuminating examples. 
A class of the so called graph states used in one-way quantum computation model [7] and fault-tolerant 
topological quantum computation 21 model as well are both good examples. Multipartite entangled pho- 
tons are essential resources used in quantum key distribution protocols |S], [2]. For the purposes of quantum 
communication a multipartite entangled states can serve as quantum communication channels in most of the 
known teleportation protocols [10]. Multipartite entanglement displays one of the most fascinating features 
of quantum theory that is called as nonlocality of the quantum world [13] . 
" Comparing to the best known two-partite case [4], [5], [6], the situation concerning a general multipartite 

, system is much more complicated even on the level of classification of the possible entanglements that can 

' arise in such systems. Given a multipartite state, a natural question to ask is whether it is entangled state 

and if so how to measure the amount of entanglement included. 

The main aim of the present note is to demonstrate how one can use exclusively two-partite systems 
methods to answer the first question for a general multipartite quantum state, see also [18) . [19] for some 
recent references on this topics. The presented in this note methods are being implemented and applied to 
study several particular (low dimensional mainly) systems within constructed in Zielona Gora University in 
the past few years simulator of quantum computation machine, called Quantum Computing System (QCS) 
which is presented partially in [20] and will be presented in [14] with all details necessary for this. The 
report from our computer experiments will be presented elsewhere. 

The present note is organized in the following way. In the next section basic definitions will be introduced 
and the algorithm detecting and localizing entanglements in the case of pure states will be formulated. An 
extension of this method to general multipartite mixed states will be presented in Section 
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2 A general multipartite systems. The case of pure states 

Let Q be a quantum system composed from n smaller and separated subsystems. The corresponding Hilbert 
spaces Hi of the subsystems are assumed to be finite dimensional with the corresponding dimensions ki for 
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i = 1, . . . ,n. 

Then the resulting Hilbert space H describing possible states of the total system is equal to H = 
"Hi 8 'H2 ® . . . "Hn and has dimension K — ki ■ k2 ■ ■ ■ ■ k„. The space of all quantum states i.e. the space of 
all nonnegative, of trace equal to one endomorphisms of the corresponding Hilbert space of states H will be 
denoted by E{'H) and its boundary dE{T-i) is composed from pure states of the system Q. 

Let p(n) stands for the set of all partitions of the set /„ = {1,2, ...,n}. For a particular tt = 
[X\,...,Xk) £ pin) we denote the number of elements of tt, |7r| = k. K partial semiorder -<j in pin) 
is introduced by the following definition: we say that tt = iX\, X2, ■ ■ ■ , Xk) is finer than a partition 
TT = (Xi, X2, ■ ■ ■ ,Xi), denoted as tt X/ tt iff for any i G {1,2,. ..,fc} there exists j £ {1,2,...,/} such that 
Xi C Xj. The maximal element vTmax exists in the poset (p(n), -<f) and TTmax = (In) and the corresponding 
minimal partitions TTmax = ({1}, {2}, . . . , {n}). 

The cardinality of the set p(n) denoted as |p(n)| is given by the so called Bell number B(n). The number 
of all partitions of the set /„ of length k is given by Stirling number of second kind S{n, k). In particular 
S{n, 2) = 2"-i - 1 and B(n) = Y.l=i Sin, k). 

With these preparatory definitions we can pass to the basic notions. 

A pure state I?/)) £ dEiHi ® "^2 (X) . . . ® Hn), n > 2 is called 7r-separable iff 

• IV') = OLilV'i} for TT = {Xi,X2,...,Xk) e p{n), It/),} G HiXi) and where 'H(Xi) stands for the 
Hilbert space corresponding to the Xi-piece of the whole system, 

• there is no finer partition tt for which \ip) is tt -separable. 

Remark 1. There is some abuse of notation introduced tn the above definition which however can be easily 
overcomed when passing to the corresponding matrix representation of 

In particular a state is called completely separable iff \ip) is TTmin-separable. Any state \tp) which is 
not TTmin-separable is partially entangled state as the following definition state. 

Definition 1. A state \tp) G dEi^ii ® 'H2 ® . . . (8 T-Ln) is caled n-entangled state iff there exists n G pin), 
TT 7^ 7i"niin such that is TV -Separable. In particular a state ji/j) which is ■Kma.^-separable will be called 
completely entangled state. 

The first question when dealing with a particular state \'4>) is the question whether this state is separable 
or entangled state. 

In the case of two-partite systems and for vector states the complete answer to this question is provided 
by the corresponding Schmidt decomposition [3], [3]. However, in general n-partite case there is no canonical 
notion of Schmidt decomposition as is well known. 

The algorithm formulated below gives a systematic way to answer this question in general and additionally 
localises the finest partition for which the separability holds. 

To start with, we introduce the following construction, which can be easily implemented in any computer 
system which contains any of the standard Linear Algebra Package. 

The canonical Schmidt decomposition: Let Hi, Hi be two Hilbert spaces of dimension n\, resp. 712 ■ 
Then for any I?/;) £H = Hi ® W2, there exists an unique decomposition: 

r(|^|>) 

J2 (1) 

i = l 

where the number 1 < ri\tp)) < min(ni,n2) is the Schmidt rank of \tp) and {\tp"),i = 1,2,3, .. . ,na} form 
complete orthonormal systems (CONS) in the corresponding spaces. 



Function SD 

input : n, ni, n2, {ip) where n — ni ■ n2, G C'"!'"^ 

output: [rAm:=U2.....,n,Am)^=1.2,...J 

remark: [r, (m) ^=u2,...,n,^ (l^')).=i,2,...,„J = ^^dV'), (m, n^)) 



We will also use the following function which computes the complete list of all 2-partitions of a given set 

X = {xi,X2, ... , X„}. 

Let us mention that it is not quite trivial to formulate the algorithm which computes the list of desirable 
two-partitions of a given n-th element set X = {xi,X2, . . . ,x„}. However using the methods presented in 
[12| this goal can be obtained. Now we are ready to formulate heuristic version of our algorithm (depicted 
as Alg. (HI)). 
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Function 2p-Par 

input : X ^ {xi,X2,---,Xn} 

output: 2pp, the list of all 2-partitions of X 

remark: [{{Xi, Zi), i = 1,2, ...)]= 2p - Par{X) 



input : = 7^1 (g) <8) • • • <8) Hn, n = k\ ■ ■ ■ ■ - kn, & dE{Hi (g) iJ2 <8> • • • <8> ^n) where fc, = dim Hi, 

output: TT G p(n) such that \tp) is 7r-separable, in particular if tt = TTmin then \tp) is completely separable. 
Algorithm 

setX = /„ = {l,2,...,n};7r=[] 
while X| > 1 do 
set N = \X\ 
for fc = 1 : AT do 

if \X[k] \ = 1 then 
X = X\X[k] 
n=[n,X[k]] 
end 
end 

set N = \X\ 
for k=l:N do 

2pL[k] = 2p-Par(X[k]) 

set i = 1 

while / < \2pL[k] \ do 
for (y, Z) G 2pL[k] do 

irY,z, |Vr>|Vz)) = SDil^Pxik]), {Y, Z)) 
if rY.z = 1 then 
set X[k] = {Y,Z) 
stop for 
end 
else 

go to the next {Y, Z) e 2pL[k] 
1 = 1 + 1 
end 

end 
end 

if I = \2pL[k] \ then 
X = X\X[k] 
7r=[n,X[k]] 
end 
end 
end 

Algorithm 1: Separable or Entangled? The case of vector states 
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3 A general multipartite system. The case of mixed states 



A general state p £ E{'Hi^y,2 ^ ■ ■ - ^T-Ln) will be called 7r-separable state, for some tt = [Xa ,Q. = 1, . . . ,k) £ 
p{n) , iff there exists a representation of p of the form 

P = y"-P« ® P^, (2) 

a 

where p^ £ EiJ-L^) for each a and moreover p is not tt -separable for any tt finer that tt. 

There exists plenty of different criterions answering (mostly partially only) the question whether a given 
state of bipartite system is separable or entangled. Without specifying the concrete form of the criterion 
used we assume that an appropriate tool for this purpose is available. It will be named as 2p-Oracle. 

2p-Oracle function 

input ■.H = Ha®'Hb where p € E{Ua ® Ub) 
output: YES, if p is entangled, NO if p is separable 



Remark 2. The best known examples of such oracles although working only for a particular classes of states 
are: PPT- criterion, certain witness construction, see f^. In general the problem of answering the 

question whether a given quantum state p is entangled or not is NP-Hard problem \17^ . 

Now, we are ready to formulate our algorithm (depicted as Alg. ([2])) by which we can answer the question 
whether a given general state p £ E{'Hi ® 'H2 ® . . . ® Hn) is completely or partially separable and in the 
second case the localization of the corresponding n G p{n) is obtained as well. 

4 Summary 

A general method for detecting and localizing partial and complete as well entanglement for general multi- 
partite states is presented. 

It was shown that the question of detection and localizing (partial) entanglement of quantum states de- 
scribing multipartite systems can be answered by using exclusively tools worked out for the two-partite case. 
However the serious drawback of methods introduced here is the non-polynomial computational complexity 
of our algorithms presented. A rough estimate gives quickly the computational complexity of both algo- 
rithms presented as 0(2"). The source of such big complexity is the length of the list of all two-partitions 
of a given set X. 

However, for small quantum registers composed of small number n e.g. 20 - 25 quantum logical units 
like qubits or qudits which is the case of all quantum computer simulators available at present (see [TJ) our 
methods might be quite useful. 

The methods introduced in the present contribution could find straightforward applications to several 
problems of quantum information processing tasks. Let us mention two of them: 

• in the process of tracing of entanglement when executing quantum programs on Quantum Computer 
Simulators, 

• in the search for a new effectively simmulable quantum circuits by using the methods of |15| . 

A natural semiorder relation with respect to the amount of entanglement included is that induced by the 
action of (S)LOCC-class of operations. In the case of two-partite systems this notion seems to be quite well 
understood [1], 0, although several basic question are still to be resolved even in this case. However, in 
the case of multipartite systems several problems arise at the very beginning. For example, let us consider 
two quantum states pi, p2 of a composite nth-partite system and let us assume that pi is 7ri-entangled and 
P2 is 7r2-entangled where of course tt; — (XJ, X2, . . . , -'ffe. ) ,i = 1,2 are the corresponding partitions of the 
set of indices {1, 2, 3, ... , n}. 

From the very definitions of (S)LOCC class of operation we must restrict ourselves to the situation 
TTi ^/ 7r2 or 7r2 ^/ tti in order to define the local actions on both states pi and p2. So, let us assume 
that TTi = {Xi,X2, . . . , Xk) -< n2 = (5^1,^2, . . . , Yi). By (S)LOCC-class of operations we mean local with 
respect to the finer decomposition {Xi, X2, ■ ■ . , Xk) completely positive operators together with general 
but local measurement equipments (supported with classical communications channels) [18) . We say that 
pi ~< p2 iff there exists (S)LOCC-class operations that transforms p2 into the state pi. From the above 

S{LOCC) 
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input 


■H 


= Ui®U2® ■ ■ - ^ Hn, n = ki-k2- ...kn, IV") G dE{Hi (g) H2 (S> ■ ■ ■ <S> Un) where h = dim Wj, 








Hi 




output 


TT e p(n) such that p is 7r-separable 


Algorithm 
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In 
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while X| > 1 do 
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1 : do 
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= 1 then 
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Algorithm 2: Separable or Entangled? General n-partite states 
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definitions it follows that we should expect much linear chains in POSET {E{Hi ® ^2 ® • . • ® Hn) , -< ) 

S(LOCC) 

then in the bipartite case. In particular it follows from [16] that there exists many locally maximal states 
with respect to this semi-order. 
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